function [xi, eta, coefs] = load_parameters(outfile)
    dresults = load(outfile, 'dresults').dresults;
    sresults = load(outfile, 'sresults').sresults;

    xi  = dresults.xi;
    eta = sresults.eta;

    coefs         = struct;
    coefs.rho     = dresults.rho;
    coefs.Pi      = dresults.Pi;
    coefs.kappa   = sresults.kappa;
    coefs.gamma   = sresults.gamma(1:sresults.ngamma);
    coefs.alpha   = dresults.theta(1);
    coefs.sigma_D = dresults.theta(2:end);
    coefs.sigma_S = sresults.gamma((sresults.ngamma+1):end);
end
